Exploring the Dark Energy Redshift Desert with the Sandage-Loeb Test 
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We study the prospects for constraining dark energy at very high redshift with the Sandage-Loeb 
(SL) test - a measurement of the evolution of cosmic redshift obtained by taking quasar spectra at 
sufficiently separated epochs. This test is unique in its coverage of the "redshift desert" , correspond- 
ing roughly to redshifts between 2 and 5, where other dark energy probes are unable to provide useful 
information about the cosmic expansion history. Extremely large telescopes planned for construc- 
tion in the near future, with ultra high resolution spectrographs (such as the proposed CODEX), 
will indeed be able to measure cosmic redshift variations of quasar Lyman-a absorption lines over 
a period as short as ten years. We find that these measurements can constrain non-standard and 
dynamical dark energy models with high significance and in a redshift range not accessible with 
future dark energy surveys. As the cosmic signal increases linearly with time, measurements made 
over several decades by a generation of patient cosmologists may provide definitive constraints on 
the expansion history in the era that follows the dark ages but precedes the time when standard 
candles and rulers come into existence. 



I. INTRODUCTION 

Measurements of luminosity distance to Type la su- 
pernovae (SNe; [1]) in combination with the location of 
the acoustic peaks in the Cosmic Microwave Background 
(CMB) power spectrum [2], as well as the scale of the 
baryon acoustic oscillations (BAO) in the matter power 
spectrum [3] provide an accurate determination of the ge- 
ometry and matter/energy content of the universe. These 
measurements are almost exclusively sensitive to the cos- 
mological parameters through a time integral of the Hub- 
ble parameter (or, the expansion rate). Although direct 
measurements of the Hubble parameter are feasible, they 
are typically difficult. For instance, the BAO probe ra- 
dial modes are directly sensitive to H(z) [4], but require 
exceedingly precise knowledge of individual galaxy red- 
shifts. In addition, the few proposals to directly measure 
the expansion rate all propose to determine H(z) at a 
few specific epochs: z < 2 (from, say, the BAO, or by 
measuring the relative ages of passively evolving galaxies 
[5]), z ~ 1000 (from the CMB [6]) and z ~ 10 9 (from the 
Big Bang Nucleosynthesis). In particular, a new cosmo- 
logical window would open if we could directly measure 
the cosmic expansion within the "redshift desert" , ideally 
exploring 2 < z < 5. 

During the early years of Big-Bang cosmology Al- 
lan Sandage studied a possibility of directly measuring 
the temporal variation of the redshift of extra-galactic 
sources [7] . As explained in the next section, this tempo- 
ral variation is directly related to the expansion rate at 
the source redshift. However, measurements performed 
at time intervals separated by less than 10 7 years would 
have failed to detect the cosmic signal with the technol- 
ogy available at that time [7]. In 1998 these ideas have 
been revisited by Loeb [8] . He argued that spectroscopic 



techniques developed for detecting the reflex motion of 
stars induced by unseen orbiting planets could be used 
to detect the redshift variation of quasar (QSO) Lyman- 
a absorption lines. A sample of a few hundred QSOs 
observed with high resolution spectroscopy with a ~ 10 
meter telescope could in fact detect the cosmological red- 
shift variation at ~ la in a few decades. In what follows 
we will therefore refer to this method as the "Sandage- 
Loeb" (SL) test. 

The astronomical community has since entertained in- 
creasingly ambitious ideas with proposals for building a 
new generation of extremely large telescopes (30 — 100 
meter diameter) [9-12]. Equipped with high resolution 
spectrographs, these powerful machines could provide 
spectacular advances in astrophysics and cosmology. The 
Cosmic Dynamics Experiment (CODEX) spectrograph 
has been recently proposed to achieve such a goal [13, 14]. 

The large number of absorption lines typical of the 
Lyman-a forest provide an ideal method for measuring 
the shift velocity. The latter can be detected by sub- 
tracting the spectral templates of a quasar taken at two 
different times. As quasar systems are now readily tar- 
geted and observable in the redshift range 2 < z < 5, we 
have a new test of the cosmic expansion history during 
the epoch just past the dark ages when the first objects 
in the universe are forming. 

If dark energy is consistent to our simplest models — 
small zero-point energy of the vacuum or a slowly rolling 
scalar field — then it significantly speeds up the expan- 
sion rate of the universe at z < 1. Under the same as- 
sumption, dark energy is subdominant at redshift z > 1, 
and almost completely negligible at z > 2. One may 
therefore ask whether is worthwhile to probe z > 2 any- 
way. The answer is affirmative: since we do not know 
much about the physical provenance of dark energy, it is 
useful to adopt an entirely empirical approach and look 
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for the signatures of dark energy at all available epochs 
regardless of the expectations. This question has recently 
been studied in some detail by Linder [15], who consid- 
ered toy models of dark energy that have non-negligible 
energy density at high redshift. 

In this paper we study the cosmological constraints 
that can be inferred from future observations of velocity 
shift and their impact on different classes of dark energy 
models. In Sec. II we review the physics behind the SL 
test, in Sees. Ill and IV we describe future constraints on 
standard and non-standard dark energy models respec- 
tively, and in Sec. V we discuss our results and future 
prospects. 



II. THE SANDAGE-LOEB TEST 

It is useful to firstly review a standard textbook cal- 
culation of the Friedmann- Robertson- Walker cosmology. 
Consider an isotropic source emitting at rest. Since it 
does not posses any peculiar motion, the comoving dis- 
tance to an observer at the origin remains fixed. Then 
waves emitted during the time interval (t s , t s + St s ) and 
detected later during (t , t Q + 5t ) satisfy the relation [7] 



dt 
W) 



t 0+ 6t a_ 
t s +st s a{t) 



(1) 



where t s is the time of emission and t Q the time of obser- 
vation. For small time intervals (St <C t) this gives the 
well known redshift relation of the radiation emitted by 
a source at t s and observed at t a , 



1 + z s (t a ) = a(t a )/a(t s ). 



(2) 



Let us now consider waves emitted after a period At s at 
t s + At s and detected later at t a + At Q . Similarly to the 
previous derivation the observed redshift of the source at 
t + At D is 

1 + z s (t D + At ) = a(t a + At a )/a(t s + At s ). (3) 
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FIG. 1: Cosmic velocity shift as function of the source redshift 
in a flat universe for different values of Qm and w. A time 
interval of 10 years has been assumed. The signal primarily 
depends on Q.m and more weakly on w. 



Friedmann equation to relate d to the matter and energy 
content of the Universe we finally obtain (see [8]): 
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(6) 



where Ho is Hubble constant, E(z) = H(z)/H is the 
(scaled) Hubble parameter at redshift z, c is the speed of 
light, and we have normalized the scale factor to a(t a ) = 1 
and neglected the contribution from relativistic compo- 
nents. For a constant dark energy equation of state we 
have 



Therefore an observer taking measurements at times t 
and t + At would measure the following variation of the 
source redshift 
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In the approximation At/t <C 1, we can expand the ratio 
a(t a + At )/a(t s + At s ) to linear order and further using 
the relation At = [a(t )/a(t s )]At s (as it can be easily 
inferred from Eq. (1)) we obtain 
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This redshift variation can be expressed as a spectro- 
scopic velocity shift, Av = cAz s /(l + z s ). Using the 



E(z) = [0 M (1 + zf + Q DE (1 + z f( 1+w ) + Q K (1 + zf 

(7) 

where and Qde are the matter and dark energy den- 
sity relative to critical, SIk = 1 — — ^de is the cur- 
vature, and w is the dark energy equation of state. 

In Figure 1 we plot Av as function of the source red- 
shift in the flat case for different values of and w 
assuming a time interval At Q — 10 years. As we can see 
Av is positive at small redshifts and becomes negative at 
z > 2. Under the assumption of flatness the amplitude 
and slope of the signal depend mainly on £1 m , while the 
dependence on w is weaker. 

In spite of the tiny amplitude of the velocity shift, the 
absorption lines in the quasar Ly-a provide a powerful 
tool to detect such a small signal. As already pointed out 
in [8] the width of these lines is of order ~ 20 km/s, with 
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metal lines even narrower. Although these are still a few 
orders of magnitude larger than the cosmic signal we seek 
to measure, each Ly-a spectrum has hundreds of lines. 
Therefore spectroscopic measurements with a resolution 
R > 20000 for a sample of ~ 100 QSOs observed - 10 
years apart can lead to a positive detection of the cosmic 
signal. Moreover astrophysical systematic effects such as 
peculiar velocities and accelerations can lead to negligible 
corrections [8]. Local accelerations may indeed be more 
important, but due to their direction dependence they 
can be determined from velocity shift measurements of 
QSOs sampled in different directions on the sky. 

In [13] the authors have performed Monte Carlo simu- 
lations of Lyman-a absorption lines to estimate the un- 
certainty on Av as measured by the CODEX spectro- 
graph. The statistical error can be parametrized in terms 
of the spectral signal-to- noise S/N, the number of Ly-a 
quasar systems Nqso: an d the quasar's redshift 
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(at z < 4) where S/N is defined per 0.0125A pixel. The 
source redshift dependence becomes flat at z > 4. The 
numerical factor slightly changes with the source redshift, 
varying from 1.4atz = 4to2atz~2. This small vari- 
ation arises because the number of observed absorption 
lines decreases with zqso- For simplicity we assume its 
value to be fixed to 1.4. The large S/N necessary to 
detect the cosmic signal implies that a positive detec- 
tion is not feasible with current telescopes. However the 
CODEX spectrograph is currently being designed to be 
installed on the ESO Extremely Large Telescope; a ~ 50 
meter giant that can reach the necessary signal-to-noise 
with just few hours integration. 

In the next sections we describe the cosmological win- 
dow that such observations could open, with particular 
focus on dark energy models. 



III. COSMOLOGICAL PARAMETERS AND 
STANDARD DARK ENERGY MODELS 

We forecast constraints on cosmological parameters 
from velocity shift measurements using the Fisher ma- 
trix method for a LCDM fiducial cosmology. We assume 
experimental configuration and uncertainties similar to 
those expected from CODEX [13, 14]. Namely we con- 
sider a survey observing a total of 240 QSOs uniformly 
distributed in 6 equally spaced redshift bins in the range 
2 ^ z < 5, with a signal-to-noise S/N — 3000, and the 
expected uncertainty as given by Eq. (8). Since there is 
no time integration involved in the computation of Av 
the Fisher matrix components can be easily determined 
analytically. 

In Figure 2 we show the near-future status of measure- 
ments of the Hubble parameter at different cosmological 
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FIG. 2: Fractional accuracy in the measurement of the Hubble 
parameter H(z) as a function of redshift for a sample of cos- 
mological probes. The accuracy and redshift ranges shown are 
best-guess values for the future surveys, and assume a single 
measurement of the Hubble parameter for each probe. Since 
SNe, number counts of clusters, and weak gravitational lens- 
ing do not measure the Hubble constant directly, but rather 
some combination of distances and (in the case of the lat- 
ter two) growth of density perturbations, we only indicate 
their redshift range with the shaded region. The Sandage- 
Loeb overall constraint assumes roughly a 30 year survey and 
other specifications as in the text. 



epochs. We assume upcoming measurements of Hq accu- 
rate to 2%, the overall BAO measurements of the expan- 
sion rate to 1% [16], the CMB measurement of 1.4% [6], 
and the SL test measurement of H using the assumptions 
outlined above. For clarity we do not show the Big Bang 
Nucleosynthesis constraint on H(z) with z = Zbbn ~ 10 9 , 
which is accurate to ~ 10% (e.g. [17]) and will improve 
as soon as the baryon density and deuterium abundance 
are more accurately determined. Since Type la super- 
novae, number counts of clusters, and weak gravitational 
lensing do not measure the Hubble constant directly, we 
only indicate their approximate redshift range with the 
shaded region. Clearly, the SL test is probing an era not 
covered with any other reliable cosmological probe. 

In Figure 3 we plot the la contours in the Qm — w plane 
as expected from type la supernovae, weak lensing power 
spectrum, power spectrum plus the bispectrum, and con- 
straints expected from Planck's measurement of distance 
to the last scattering surface. Supernova and weak lens- 
ing estimates are based on the SNAP mission [18], and 
SNe include systematic errors. The SL contour are com- 
plementary to those of other tests since Av probes a dif- 
ferent degeneracy line in the Qm ~ w plane. It is worth 
noticing that such measurements are mostly sensitive to 
the matter density, as expected from the plot shown in 
Fig. 1. Allowing also for variation of the curvature, fix, 
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FIG. 3: Constraints in the Qm — w plane for the SL test 
assuming (with increasingly smaller contours) a 10, 30 and 
50-year survey and the number of quasars as in the text. We 
also show constraints expected from type la supernovae, weak 
lensing power spectrum, power spectrum in combination with 
bispectrum, and constraints expected from Planck's measure- 
ment of distance to the last scattering surface. SNa and weak 
lensing estimates are based on the SNAP mission, and SNe 
also include systematic errors. The SL test contour is comple- 
mentary to other constraints, and is mostly sensitive to the 
matter density, represented here by Qm- 



> 

< 




FIG. 4: Cosmic velocity shift as function of the source redshift 
for a fiat ACDM model with Q m = 0.3 (solid line), Chaplygin 
gas model (long dashed line) and an interactive dark energy- 
dark matter model (short dashed line) with 2% deviation from 
standard redshift scaling of matter (see text). The error bars 
on the ACDM model assume a 10 year survey and other spec- 
ifications as in the text. 



IV. NON-STANDARD DARK ENERGY 
MODELS 



we find that the SL test alone determines the matter den- 
sity about four times better than the curvature. As an 
example, for the 30-year survey the marginalized errors 
are ct(Q,m) — 0.03 and a (fix) = 0.13. 

We have also found that limited accuracy in the Hubble 
constant measurement does not degrade the power of the 
SL test. For example, if h is known to 0.04 (or to about 
5%), the accuracy in w is degraded by only 2% relative 
to the case when h is perfectly known. 

As far as dark energy is concerned, assuming a Gaus- 
sian prior on h with ah = 0.04, we obtain a w — 0.8 for 
Afo = 10 years time interval and 0.3 for 30 years. Thus 
the constraints on w are not competitive with those in- 
ferred from other, well-established tests such as SNe la, 
weak lensing or BAO once we take into account that the 
latter probes will provide strong constraints by the time 
the SL test is undertaken. However, one should note 
that the constraints obtained by SL decrease linearly with 
time. For measurements made over a century, and with 
the expectedly larger number of QSOs, the SL limits on 
w can easily be at the few percent level. 



A unique advantage of the SL test is that it probes 
the redshift range 2 < z < 5 which is very difficult to 
access otherwise. As mentioned in Section I, probing this 
redshift range is important for testing non-standard dark 
energy models that would otherwise be indistinguishable 
from those with a smooth, nearly or exactly constant 
equation of state function w(z). 

WMAP observations in combination with low-redshift 
limits from SNe la impose a weak upper bound on the 
amount of dark energy deep during matter domination 
to be Qde(z) < 0.1 [19, 20]. Let us suppose that dark 
energy re-emerged at 3 < z < 5 with Ode(^) = 0.1 
in this range, while essentially zero at larger redshifts 
and that it behaved as a standard ACDM at z < 3 (i.e. 
assuming standard Ode = 0.75, w — —1 values). We 
would like to distinguish such a model from a pure ACDM 
with Ode and w as above. Clearly, low-redshift probes 
(SNe, BAO, galaxy clustering) cannot distinguish these 
two models since the required redshift is too high even for 
the most ambitious surveys. Furthermore, the distance to 
the last scattering surface between the two models differs 
only by 0.5%, which is too low to be observable even 
with Planck, since the 1-a uncertainty is about 0.4% with 
temperature and polarization information [21]. However, 
the two models can be distinguished via the SL test at 
about 3-cr level, assuming only a 10-year survey and other 
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specifications as in the previous section. 

While the aforementioned scenario with dark energy 
emerging in the specific window at 3 < z < 5 may seem 
contrived, it is easy to find physically motivated models 
whose identification can significantly benefit from data in 
this "desert" . One example is given by scalar field mod- 
els which predict the periodic emergence of DE at var- 
ious epochs during the history of the universe [22, 23]. 
Even though for these particular models one has to go 
to a much higher redshift to see the next phase where 
^de = 0(1), it is plausible, and certainly currently ob- 
servationally allowed (e.g. [24]) that such a phase could 
have occurred somewhere within 2 < z < 5. 

Another example is given by models with dark energy- 
dark matter interaction (see e.g. [25-31]). In the simplest 
realization where the scalar field only couples to dark 
matter, it mediates a long range interaction which causes 
two separate effects. First, dark matter particles, unlike 
the baryons, experience a scalar-tensor type of gravity, 
which modifies the Newtonian regime (see [32]). The 
time and scale of when such type of modification becomes 
cosmologically relevant depend on the particular model 
considered. Second, dark matter particles acquire a time 
dependent mass whose evolution is determined by the 
specifics of the scalar field dynamics. As a consequence 
of this, the redshift evolution of the dark matter density 
deviates from the usual (l+z) 3 . At low redshift, when the 
universe is dark energy dominated, these models cannot 
be distinguished from the standard ACDM . Therefore, 
the best way to probe models with such dark matter-dark 
energy interaction is to map out cosmic expansion during 
the matter dominated phase (see Figure 4). The SL tests 
offers a unique tool to do just that. 

In order to forecast how well a deviation of the dark 
matter density from the (1 + z) 3 law can be detected, we 
parametrize its redshift evolution as oc (l + zf^-V in the 
range 2 < z < 5, where |6| < 1 is a constant free param- 
eter. The scalar field, on the other hand, can be treated 
as a dark energy component with w w — 1, since the field 
slowly rolls toward the minimum of its effective poten- 
tial at late times [30]. Assuming a flat fiducial model 
with Um — 0.3 and w = —1, we find that the SL test 
can detect deviations from the standard matter scaling 
as small as 1% (i.e. b = 0.01) over 10 years and 0.3% for 
30 years. Therefore, SL test can provide constraints an 
order of magnitude tighter than those inferred using fu- 
ture SNe la or the Alcock-Paczynsky test [33] . Since the 
deviation b is generally a function of redshift, one can use 
the velocity shift measurements to reconstruct the red- 
shift dependence of b, and then determine the strength 
and functional form of the scalar interaction. 

The Chaplygin gas is yet another dark energy candi- 
date that can be tested in the range of redshift probed by 
the SL test. Proposed as a phenomenological prototype 
of unified dark energy and dark matter model [34-36] , it 
describes an exotic fluid with an inverse power law ho- 
mogeneous equation of state, P = — \wo\il^, h pQ / p a (e.g. 
[37]), where wq is the present equation of state, fip h is the 



current energy density of the gas pa = 3Hq is the total 
energy density and a > —1 is a dimcnsionless parame- 
ter. This corresponds to a fluid which behaves as dust in 
the past and as cosmological constant in the future. For 
a = the model reduces to ACDM [38]. The Chaplygin 
gas energy density evolves with redshift according to: 



n Ch (z) = n° Ch [-two + (i + «*>)(i + ^) 3(a+1) 



(9) 



As shown in [39] this model can provide a good fit to 
current cosmological observables with fich ~ 0.95 (with a 
baryonic component of Ot. = 0.05), wo ~ —0.75 and a = 
0.2. From Fig. 4 we can see that although this model has 
a velocity shift at z < 1 similar to that of a ACDM , it can 
be tested with high redshift measurements. For instance 
we find that, for this particular model, the Chaplygin 
parameters can be determined with uncertainties a Wo = 
0.03 and a a = 0.04 respectively and thus distinguished 
from the ACDM values at a high confidence level. 



DISCUSSION 



In this paper we have analyzed the prospects for con- 



straining dark energy at high redshift (2 
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by direct measurements of the temporal shift of the 
quasar Lyman-a absorption lines (the Sandage-Loeb ef- 
fect). While the signal is extremely small, the physics is 
straightforward, and the measurement is certainly within 
reach of future large telescopes with high resolution spec- 
trographs. 

As the SL test mostly probes the matter density at high 
redshift, the constraints on standard dark energy mod- 
els with a nearly or exactly constant equation of state w 
are weaker than those that observations of SN la, BAO, 
weak lensing and number counts will be able to achieve 
in the future (although the SL test becomes competi- 
tive for measurements spread over a period of several 
decades or more). This is mostly because the sensitivity 
of cosmological probes to standard dark energy models 
is exhausted at z < 2 and higher redshift data do not 
improve them significantly (e.g. [40]). However, the prin- 
cipal power of SL measurements at 2 < z < 5 comes from 
their ability to constrain non-standard dark energy mod- 
els, where the dark energy density is non-negligible at 
higher rcdshifts or equivalcntly, models where the to- 
tal energy density does not scale with redshift as (1 + z) 3 
at z > 2. In particular, as discussed in the previous sec- 
tion, in only one decade the SL measurements will allow 
to test the redshift scaling of the matter density with an 
accuracy one order of magnitude greater than standard 
cosmological tests. 

Tighter constraints on standard dark energy models 
could be obtained if observations of Ly-a systems were 
feasible at z < 2. What are the prospects for perform- 
ing the SL test at lower redshift? For this to be possible 
UV space-based instruments are necessary. Space-based 
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UV Lyman-a astronomy is possible and has already pro- 
duced remarkable results (see e.g [41]), though it gener- 
ally lacks the spectral resolution and wavelength coverage 
of the higher rcdshift studies. However during the past 
few years high quality observations of several low red- 
shift QSOs have been obtained with the Hubble Space 
Telescope (HST) and its Space Telescope Imaging Spec- 
trograph (STIS). It is therefore conceivable that future 
space based experiments will be able to measure the SL 
effect at low rcdshift. 

Other astrophysical probes that can potentially ex- 
plore the redshift desert are not yet well understood. For 
instance, it might be possible to measure the angular di- 
ameter distance at redshift of 10-20 from the acoustic os- 
cillations in the power-spectrum of the 21cm brightness 
fluctuations [42]. Further, gamma ray bursts (GRBs) 
(e.g. [43]) have been proposed as alternative standard 
candles. These can probe roughly the same redshift range 
as the SL (z < 6). Similarly gravitational wave "standard 
sirens" can measure distances to GW sources all the way 
out to z <~ 20 [44]. However, the prospects of turning 
the GRBs into standard candles is extremely uncertain 
at this time, and it is not clear that they can be used as 
cosmological probes. The GW sirens, while potentially 
providing amazingly accurate distance measurement, suf- 
fer from gravitational lcnsing of the signal, which adds an 
effective noise term to distance measurements and greatly 



degrades their power to constrain cosmological models 
[44]. Conversely, the SL test is based on extremely sim- 
ple physics and involves controllable systematic errors, 
but it docs require a powerful instrument and patience 
to wait at least a decade before repeating the measure- 
ments in order to produce interesting results. 

We finally point out that the velocity shift sig- 
nal increases linearly with time, thus amply rewarding 
increased temporal separation between measurements. 
Therefore observations made over a period of several 
decades by a generation of patient cosmologists may pro- 
vide definitive constraints on the expansion history in the 
era before the usual standard candles and rulers, Type la 
supernovae and acoustic oscillations in the distribution of 
galaxies, become readily available. 
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